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Abstract 

In [2], F.Linares, G. Ponce, J-C.Saut have proved that a non-fuUy 
dispersive Zakharov system arising in the study of Laser-plasma inter- 
action, is locaUy well posed in the whole space, for fields vanishing at 
infinity. Here we show that in the periodic case, seen as a model for 
fields non- vanishing at infinity, the system develops strong instabilities 
of Hadamard's type, implying that the Cauchy problem is strongly 
ill-posed. 
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1 Introduction 

The construction of powerful lasers allows new experiments where hot plasma 
are created in which laser beams can propagate. The main goal is to sim- 
ulate in a laboratory nuclear fusion by inertial confinement. We need some 
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precise and reliable models of laser-plasma interactions in order to produce 
numerical simulations than are useful in order to understand the experi- 
ments. The kinetic-type models are the more precise ones but the cost in 
term of computations is exorbitant and no physically relevant situation for 
nuclear fusion can be simulated in this context. Therefore, we use a bi- 
fluid model for the plasma : we couple two compressible Euler systems with 
Maxwell equations. Even under this form, it is not possible to perform direct 
computations because of high frequency motions and the small wavelength 
involved in the problem. At the beginning of the 70's, Zakharov and its 
collaborators introduced the so-called Zakharov's equation [7| in order to 
describe the electronic plasma waves. Basically, the slowly varying envelope 
of the electric field E = ^tp is coupled to the low-frequency variation of 
the density of the ions 5n by the following equations which are written in a 
dimensionless form : 



(1.1) 



' ^9*^ + A(VV') = VA-idiv(5nV^A), 
_ 5|(5n - A(5n = A(|V^|2). 



Of course, variations of this systems exists (see [6] for example). For laser 
propagation, one uses the paraxial approximation and the Zakharov system 
reads 



(1.2) 



i{dt + d,)E + A^E = nE, 
{d^ - A,)n = A,\E\^ 



The space variables are {z,x), z S M and x £ M?; z is the direction of 
propagation of the laser beam and x are the directions transversal to the 
propagation. See [5] or j6] for a symmetric use of this kind of model for 
numerical simulation. 

We consider the Cauchy problem for ()1.2p with initial data 



(1.3) 



E\t=o — Eq, 

n\t=o = no, dtn\t=o = rii. 



The existence theorem (see [6l IT] [3j and references therein) for the classical 
Zakharov system, that is when A^; is replaced by A(j,,a;)) does not apply. In 
[2], it is proved that the Cauchy problem for (jl.2p is well posed, locally in 
time, for data in suitable Sobolev spaces. The proof is based on dispersion 
estimates. For periodic data, these dispersion estimates are not valid. This is 
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a well known phenomena, even in the simpler case of Schrodinger equations. 
However, the new phenomena here is that the consequences of this lack of 
dispersive effects are much more dramatic since it implies strong instabilities 
of Hadamard's type, so that the Cauchy problem for periodic data is strongly 
ill-posed in Sobolev spaces. 

For numerical experiments, one uses mainly periodic boundary condi- 
tions: one considers that the plasma is infinite and has a periodic structure. 
For this kind of application, it is quite reasonable to consider that E ^ 
at infinity. Our result has therefore a practical application and means that 
the paraxial approximation is not a good model in this case: one should add 
the longitudinal dispersion. 

We look for solutions U = {E,n) of (jl.2p . which are periodic in x, with 
period 27r in x and periodic in z with period 2'kZ, where Z is arbitrary. We 
denote by T the corresponding torus M/27rZ x (]R/27r)^. 

We consider the constant solution 

(1.4) U={E,Q), E^O, 

which of course does not belong to the spaces used in [2] , and we prove that 
this solution is strongly unstable. 

Theorem 1.1. For all s, there are families of solutions Uk = U_ + {&k,nk), 
in C^{[0,Tk];H''{T)) such that 

(1.5) ||efc(0), nfc(O), 5*71^(0) ||j^.(t) ^ 0, 

(1.6) Tk ^ 0, 

(1-7) ||efc(Tfc), nk{Tk)\\L2(r) ^ oo. 

This nonlinear instability result is pretty strong: not only the amplifica- 
tion \\u{T)\\o/\\u(0)\\s is arbitrarily large, in arbitrarily small time T, with 
arbitrary loss of derivatives s, but there is an effective blow up of the 
norm. 



2 The instability mechanism 

Our construction is based on the analysis of the dispersion relation for the 
Zakharov system. Consider the linearized equations around {E, 0) : 




i{dt + dz)e + A^e - En = f, 
{df-A^)n-A^{2ReEe)=g 
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With (e,e,n) as unknowns the system reads: 

-i{dt + d,)e- A^e + En = -f, 

(2.2) I i{dt + d^)e - Ae + En = -J , 

, {d'f - A^)n - ^ A^e - 1 A^e = g. 
Denoting by (t, C,0 the frequency variables dual to (t, z, x), its symbol is 

/(r + c) + iep 

(2.3) 0_ -(r + C) + |?P 

and the relation dispersion is P = 0, where P is the determinant of the 
system, that is 

(2.4) P = (|e|2 - r2)(|er - (r + C?) - 2|^Piei' = ^'o - 2|^|2|e|^ 
The remark is that for (C, ^) real, Pq has 4 real roots in r 

(2.5) -lei, -c-iep, -c + iep, 

with an intermediate double root when < |^| = — C — Note that Pq is 
of degree 6 in while the perturbation — is of degree 4 and negative. 
Therefore, for ^ large and Q = — 1^| — the double root of Pq is perturbed 
in two conjugated complex roots. More precisely, for 

iei»o, c = -i^i-iei' and T = iei(i + (T), 

the determinant P is 

(2.6) P = -lel' (^'(2 - a/|e|)(2 + a) + 2\Efm) . 

The implicit function theorem shows that there are two non-real roots 

(2.7) ^ = ^±i^|^||+o(l). 

This means that waves at frequency (CiO with C, = — 1^| — are am- 
plified by the exponential factor 

This implies that the Cauchy problem for the linearized equations (|2.ip is 
ill-posed in : there are Cauchy data in such that the homogeneous 
problem with f = g = has no solution in C'^([0, T]; H^°°). 

The goal of this paper is to translate this spectral instability into a non- 
linear instability result for the Zakharov system (II. 2p . 



(2.8) e^*l«l^, 7 = ^>0. 
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Remark 2.1. How is it that this spectral instabihty does not intervene in 
the analysis of [2]? The first answer is that the condition ^ 7^ is cru- 
cial for 7 to be positive. In their case, where solutions vanish at infinity, 
linearizing the equation around non-vanishing constants has no real signif- 
icance. However, the symbolic calculus above also makes sense in the case 
of variable coefficients and one expects that the dispersion relation P = 0, 
with E replaced by E{t, z, x), which still has non-real roots, should play an 
important role in the analysis. For instance, the symbolic analysis appears 
when one replaces the plane wave analysis used for constant coefficients, by 
geometric optics expansions associated to localized wave packets. In this 
case, for a wave packet with mean frequency (— 1^| — an exponential 

amplification similar to ()2.8p is expected. But the group velocity in x of 
this packet is of order 2^; therefore if E is confined (think of it as compactly 
supported) the time of amplification is short (typically 0(|^|~^) ) so that the 
overall effect of the amplification is bounded. Of course, this is just a very 
rough explanation, but it is rather intuitive. The detailed balance between 
amplification and localization is indeed given by the dispersive estimates 
proved in [2]. 

Remark 2.2. The system can be reduced to first order in t, introducing 
{dxC, dtn, dxTi) as unknowns, but it is not first order in x, because of the 
Schrodinger part of the system. However, there is a good analogy with the 
analysis of weakly hyperbolic system. Indeed, the analysis of the symbol 
(j2.3p shows that when for = — — |^|, there is a double eigenvalue with a 
2x2 Jordan block. The exsitence of non-real eigenvalues (|2.7p . simply means 
that the natural analogue of the Levi condition for first order system is not 
satisfied. Pursuing the analogy, the exponential growth ()4.4p indicates that 
the Cauchy problem should be well posed in Gevrey classes for s < 2. 



3 Scheme of the proof 

It is certainly sufficient to prove the theorem with functions of x = (xi,X2) 
independent of X2- To simplify notations, we assume from now on that x is 
one real variable. Consider spatially periodic solutions of (II. 2p . with period 
27r in x and 27rZ in z. Moreover, we look for solutions n and E of the form 

n = n(kx — mz, t) 
(3.1) ^ ^ ^ ^ 

E = E + e{kx — mz, t) 
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with new functions n{9,t) and e{6,t) 2tt periodic in 9. For the functions to 
be 2tt periodic in x and 2tt/Z periodic in z, it is sufficient that 

(3.2) ken, mZ en. 

To be close to the unstable frequencies, we require that \m — k — k'^\ <C Vk 
and therefore we choose m G n/Z such that 

(3.3) {k'^ + k)-l/Z <m<{k'^ + k). 
The new equations read 

{i{dt — md0)e + k'^dgc — En = ne, 
{dl - k^dl)n - k^dl{Ee + Ee) = k'dj\e\\ 

With U := *(e,n), write it as 

(3.5) Lk{dt,de)U = Nk{U) 

where Lk is the linear operator defined in the left hand side of (13. 4p . and 
Nk{u) the quadratic term in the right hand side. 

The first step concerns the homogeneous equation 

(3.6) LkU = 0, 

which is studied using Fourier series expansions in 6. The choice (13. 3p and 
the spectral analysis of Section 2 and the choice (j3.3p imply that for k large, 
the harmonic 1 is unstable : 

Proposition 3.1. There is ko such that for k > ko, there are solutions 

jja ^ (ga^^a^) gf guch that 



(3.7) 

with 



e" = e1{t)e'^ + eli{t)e~'^ 
= sinh [ta] cos (tRe X + 9) 



(3.8) 4i(t) = (e"±i,+e*^ + e"±i,_e-*^) 



where the parameters A, a, e±i^± depend on k, A and a being real positive 
and, as k ^ +oo, there holds : 

(3.9) e^i + ~ -lEj^a, e\^ _ ~ -iEj^o, e^^ ± = 0{k~^). 



(3.10) A~A;, cj~|^|Vfc/2. 

The proof is given in Section 4. 
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Next, we consider SU"" as a first approximation of tlie solution of (|3.5p 
to construct, with 5 a small parameter to be chosen. More precisely look for 
solutions of (13.51) as 



(3.11) U = 5{U'' + u), u={e,n) 

with the same initial data as SU"". Because the nonlinearity is exactly 
quadratic, the equation for u reads 

(3.12) Lk{du de)u = dNkiW" + u) , 6,^=0 = n\t=o = dtn\t=o = 0. 

This equation is solved by Picard's iteration and the main step is to solve 
the linear equation 

(3.13) LkU = F, e|t=o = "-14=0 = 5*^^14=0 = 0. 

in Banach spaces which are also well adapted to the nonlinearity. The choice 
of these spaces, more precisely of their norm, is technical and dictated by 
the computations detailed in the next sections. We just give here their 
definition. 

For a periodic function t> of ^, we denote by Vp its Fourier coefficients so 
that 



(3.14) v = Y.^P 



The first Fourier coefficient ei plays a special role and we use the notations 

(3.15) e{t,9) = ei{t)e'^ + e'{t,e). 

For s > 1 and T > , se denote by E.^{T) the space oi u = (e,n) with n 
real valued, such that 

(3.16) eeC^{[0,T];H'+^)nC\[0,T];H'), n e C\[0,T]; H') 
equipped with the norm 

II'"IIei(T) = sup e'^'^lk^leiit)] + k~^\dtei{t)\ + k^\e{t)\\Hs+2 

(3.17) *^[°'^] 

+ k-^\\dte'{t)\\H^ + \\n{t)\\H^ + k-^dtn{t)\\Hs] 

where a is defined at Proposition 13.11 The norm depends on k > 1 and s, 
but, to lighten the text, we do not mention this dependence explicitly in the 
notations. 
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(3.18) 



We denote by E^{T) the same space ()3.16p . equipped with the norm 

n||E2(T) = sup e-^'''{k\ei(t)\ + \dtei{t)\ + k\\e'{t)\\Hs+2 

iG[0,T] 

+k--i\\dte'{t)\\Hs +k^\\n{t)\\H^ +k~^\\dtn{t)\\H^y 



There are two differences between p.l7p and (|3.18p : first the weight e"'^* 
is replaced by e"^*^* and second all the powers of k in the coefficients are 
increased, at least by a factor i. In particular, 

(3.19) IhllEHT) < ^"^e'^'^ll^^llE2(T) 

For the right hand sides, we denote by F^(T) the space of = {f,g) 
with g real valued such that 

(3.20) f eC\[0,T];H'), g e C\[0,T]; H') with 50 = 0, 
equipped with the norm 

||F||ei(t)= sup e-"''\k"2\\f{t)\\Hs+k-^\dtfmHs 

(3.21) *^[°'^] ^ 

+ k-^\g{t)\\Hs} 

The next three results justify the choices of these norms. We assume 
that the parameter s > 1 is fixed. 

The first estimate is an immediate consequence of Proposition 13.11 and 
([3:91) (IXTnl) . 



Lemma 3.2. There is a constant K"^ such that for all k > ko and all T < 1, 
the approximate solution U"" of Proposition \3.1\ satisfies 

(3.22) I|C^"IIei{T) < K''. 

The next two propositions are proved in Section 6. 

Proposition 3.3. There is Ci > 0, such that for all k > ko, all T < 1 and 

all F G ¥^{T), the Cauchy problem ([3T3|) has a unique solution U £ E^{T) 
and 

(3.23) l|f^llE2(T) < Ci||-^IIe2(t)- 

The nonlinearity Nii.{U) occurring in (j3.5p is quadratic. Denote by 
■^k{U, V) the bilinear associated form such that Nk{U) = Mk{U, U). 
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Proposition 3.4. There is C2 > 0, such that for all k > ko, allT < 1 and 
all U and V in E^{T), there holds Mk{U,V) G ¥'^{T) and 

(3.24) \\Mk{U,V)h2(T) < C2||C/||ei{t)II^IIei(t)- 

These estimates easily imply the following: 

Corollary 3.5. There are cq > 0, C and ko, such that for all k > ko and 
all 5 g]0, 1], the problem (j3.12p has a unique solution u = (e,n) in the unit 
ball ofK^{T), provided that 

(3.25) dk'h''^ < Co. 
Moreover, the solution satisfies 

(3.26) \Ht)\\Hs <Ck-h''^ 

Proof. Denote by L~^^F the solution of (j3.13p . and consider the mapping 

u ^ Tu:= 5L^^Nk{u'' + u) 

which, by the lemma and propositions above, is well defined from E^(T) to 
Ei(r). Moreover, 

||Tn||iEi(r) < CiC2<5A:-3e'^^(K" + ||^x||e(t))'- 

Thus it maps the unit ball to of E^(T) to itself, if (|3.25|) holds with cq small 
enough. Similarly, decreasing cq if necessary, one shows that this mapping 
is contractive on the unit ball, implying the existence and uniqueness of the 
solution of tt = Tu in the unit ball. 

The equation u = Tu and the estimates also imply that 

\\n{t)\\Hs < A:-^e2'"*||u||E2(T) < CiC25k-^2e^''\K'' + 1)^ 

< CiC2CQk~-^e''\K'' + lf 

finishing the proof of the Corollary. □ 

We end this section by proving that the main Theorem 11.11 is a conse- 
quence of this analysis. 
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Proof of Theorem \l.l\ We fix an integer s. With 
(3.27) 5 = A;-(2^+2)^ 

Corollary 13.51 provides us with solutions of (|3.5p . Uk = IJ_ + 5{U°' + Uk)-, with 
Uk in the unit ball of W}{Tk) and = ^ ln{k'^'''^'^^i /cq) satisfies 



(3.28) 6k-^/^e''^^ 



Co. 



Since a is of order by ()5.7p . T^, tends to as A; tends to infinity. 

Going back to the (z, x) variables, according to the change of variables 
()3.ip . we obtain solutions, denoted by = U_ + Uk, of the original Zakharov 
system ()1.2p . Set = (cfci^fc); these functions are deduced from 5{U^ + Uk) 
by the change of variables (j3.ip . Since m <k'^ + k, we can evaluate the ff** 
norm (in the variables {z,x)) of the Cauchy data 

\\{ek\t=o,nk\t=o,dtnk\t=o\\H'>iT) < C 5 k^^^'^ WU" + Ukiy (t) 

< CSk^'+\K'' + 1). 

Note that there is no Jacobian factor because the norms are taken for 
(z, x) G T in the left hand side and for 6 £ M/27rZ in the right hand side so 
that 

f measT f'^^ 

(3.29) / v{kx - mz)dzdx = / v{e)de. 

Jt 27r Jq 

Therefore, with our choice of S, the left hand side tends to zero as k tends 
to infinity. 

Finally we compute the norm of hk at time T^.. Using (|3.29p and 
p.7p we see that 

\\nkiTk)\\L2(j) > ci6smh{Tka) - 6\\nk(Tk)\\L2 

with ci > independent of k. Therefore, (j3.26p (j3.28p imply that 

||nfc(Tfc)||i2(Tr) > ^CK^e'^^'^ - Cdk'U'^^'' - 0{6e-''^'^) 
> ^cicok^ - Cco - o(l). 

Therefore this norm tends to +cxd and the proof of the theorem is com- 
plete. □ 
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4 The linear instability 

We study the linear equation for U = (e, n) and F = (/, g) 

(4.1) LkU = F 
using Fourier series expansions in 9: 

(4.2) e{e,t) = Y,epit)e'P', n(0, t) = np(t)e*P^ 
Since n and g are real, 

(4.3) n_p = fip, g^p = gp, 
and ()4.ip reduces to 



(4.4) 

and for p > 1 

(4.5) 

with 
(4.6) 



Fn 



' {idt + mp- k'^p^)ep - Ehp = fp, 
{idt + mp + k^p'^)ep + Ehp = fp, 
I (9| + /cV)np + k^P^CRep + Eep) = gp, 



6— P) fp — /- 



are the Fourier coefficients of e and — / respectively. For p > 0, we denote 
by L]^{dt,p) the linear operator in the left hand side of (|4.5p . 



In the remaining part of this section we concentrate on the case p = 1 
and prove Proposition 13. 1[ We reduce ()4.5p for p = 1 to a first order system 
by introducing vi = —ik~^dthi. The equation reads 



(4.7) 
with 

and 



idtVi + AVi = Fi, 
Vi = (ei,ei,ni,wi), Fi = {fi, fi,0, k'^ gi] 



A 



(m-k'^ -E 0\ 

m+fc2 ^ 

A; 

\ k'E kE k 
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Lemma 4.1. If and k is large enough, A has four distinct eigenval- 
ues; two, called Ai and A2 are real and the other two, A3 and X4, are non-real 
and complex conjugated. There holds 



(4.8) 



Ai ~ 2r 



As 



-k, ReXsr^k, a := Im A3 ~ |£:| ^^72. 



Proof. This follows from the analysis of the determinant equation in Sec- 
tion 2. The eigenvalue equation is 



P = (A2 - k^){{X - mf - k^) - 2\E\^k^ 







Following (j3.3p . we write m = A;^ + A; + m', and the equation reads 

(A^ -k^){X-k + m'){X - 2k^ - k + m') = 2\E\'^k^ 

Because m' = 0(1), the lemma easily follows by perturbation analysis of 
the roots of 

(A^ -k'^){X-k + m')(A - 2k^ - k + m') = 0. 

□ 

Next, to evaluate e**^, we need to analyze the eigenprojectors of A. 
Denote by rj [resp. Ij] right [resp. left ] eigenvectors of A associated to the 
eigenvalue Xj. Then 



(4.9) 



JtA 



AtXn 



A detailed inspection of the eigenvector equations implies the following 



n 



(o{k-^)\ 

1 

0{k-^ 

\o{k-')J 



, h = {0{k-^), 1, 0(A:-2), Oik-^)) 



r2 



/0(fc-i)\ 
0{k-^) 
1 

V -1 / 



Z2~(0(l), Oik-^), 1, -1), 
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where for vectors, a ^ h means that aU the components satisfy au ~ bk- 
Moreover, 



/ iE/a \ 



h^{kE/ia, 0{k-^), 1, 1) 



1 



V 1 / 



r4 = r3, = I3 



where = k\E\'^/2 w k. Note that r^^A = 0(1) and - u = 0{\E\/Vk) 
and ^3^4 = 0{\E\\/%) while r3^4 • ^3^4 ~ 4. This reflects that for ^ = 0, the 
corresponding matrix has a Jordan block. 

Proof of Proposition \3.1[ 
With notations as above. 



is a solution of (|4.7p with Fi = 0. It corresponds to a solution (e", e", n") of 
LiUf = and therefore to a solution 



of LiU" = 0. 

Choosing, as we may, r^, and r4 such that the third component is exactly 
equal to one, we obtain that 

n"-{t, 9) = sinh {ta) cos (tRe A3 + d) 

and the estimate ()3.9p follows from the estimates of the eigenvectors above. 
Moreover, (IS.lOp follows from Lemma l4.ll □ 

Next we turn to the analysis of (j4.7p . The solution with vanishing initial 
data is 





rf = fife'" + nfe- 



(4.10) 
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Introduce = Ij • Fi. With / denoting (/i, /i) and g = gi there holds 

^1 = *f + *k-'^g, 
(4.11) = */ + *A:-i5, 

$3 4 = *y/kf + *k-^g. 

where * denotes constants coefficients that are uniformly bounded in k. Let 



The properties of the rj 's and (|4.10p imply that the components (ei, ei, ni, fi] 
of Vi satisfy: 

(4.12) 

ni = *k + *^'2 + *^3,4, 

Vi = *k^^^i + + *^3,4- 

We use the following elementary estimates: 

Lemma 4.2. Let ^ 

^{t) = [ e^^(*-"Vj(s)ds. 

^0 

There holds 

t 



m)\ < / e-i-^(*-^)|,/.(s)|ds, 
Jo 

\dtm\<\>^j\\m\ + \mi 
\Xj\\m\ < i5tV'(t)i + i</'(t)i. 







To simplify notations, we note ^4 < S to mean that there is a constant C 
independent of k such that A < CB. We use the first and second estimate 
of Lemma 14.21 to bound the contributions of g to the integrals in (j4.1(jp . 
and we use the third and fourth estimate, when necessary, to bound the 
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contributions of /. Therefore, 

l^iWI < f\f{s),k-S{s)\ds, 
Jo 

mm < 1/(0)1 + \k-^9{t)\ + r \dtf{s), k-^g{s)\ds, 

Jo 

k'mt)\ < 1/(0)1 + |/(t)|+ f\dtf{s),k-^g{s)\ds 

Jo 

l*2(t)| < ['\f{s),k-'g{s)\ds, 
Jo 

m2{t)\ < 1/(0)1 + \k-'g{t)\ + r \dtf{s),g{s)\ds, 

Jo 

k\Mt)\ < 1/(0)1 + \f{t)\ + f \dtf{s),g{s)\ds, 

Jo 

l^3,4(t)| < f e^'-''>n^f{s),k~'g{s)\ds, 
Jo 

m3A{t)\ < e"^\Vkfm + \k-'g{t)\ + f e^'-')"\^dtf{s),g{s)\ds, 

Jo 

k\^ZA{t)\ < e'^\Vkfm + \^f{t)\ + f e^'-^^"\^dtf{s),g{s)\ds. 



Adding up the various estimates, we obtain: 

Proposition 4.3. For p = 1, the solution (ei,ei,ni) of (|4.5p with vanishing 
initial data satisfies: 

\eiit)\ < [\^^'-'^\fi{s),k-lg,{s)\ds, 







(4.13) \dMt)\<e''yim + \k-hi{t)\ 

ft 

e 





'''^'-'^\dth{s),k~hi{s)\ds, 



(4.14) 



k'\ei{t)\ + \dtei{t)\ <e'^*|/i(0)| + |/i(t)| + \k~^gi{t)\ 

+ / e''^'-'^\dtfiis),k-'gi{s)\ds 







k\Mt)\ + \dtMt)\ <e'^'k-2\M0)\ + \k-2Mt)\ + \k-'gi{t)\ 

(4.15) rt . 

+ / e^^'-'^\h{s),k-2dtMs),gi{s)\ds, 
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where fi = (/i,/i). 

Corollary 4.4. There are ko and C such that for all k > ko, K , T > 0, 

and all fi = (A, fi), gi satisfying for t £ [0, T] 

k'^\fi{t)\ + k-'^\dtfi{t)\ + k-^g,{t)\ < Ke"^\ 

then the solution of ()4.5p for p = 1 with vanishing initial data satisfies 

k\ei{t)\ + \dMt)\ <CKe^^\ 
k\e_i{t)\+k-^\dte-iit)\ < CKe^''\ 
kh\hx{t)\ + k-\\dthx{t)\ < CKe'^'^K 

Proof, a) From Proposition 14.31 we deduce that 

k\ei{t)\< CKi^ f e^^^'^'^e^'"'' dt' < CKe^"', 
Jo 

where we have used that a ^ Vk. Similarly, 

\dtei{t)\ < Ci^i(A:-V2e-* + A^-^/^gS-* + ^e'^^'-^'h^'"'' dt') < CKe^'^K 

This implies the first estimate. 

b) Similarly, ()4.14p implies that 



(4.16) 



fc2|e_i(t)| + \dte.i{t)\ < C7Ki (^e"* + e^'^* + j' y/ke"^''''^ e^"^' dt' 

< CKe^"K 



c) The estimate ()4.15p implies that 

A;|ni(t)| + \dtni{t)\ < CKiie"^ + e^-^* + f ke'^^^'^'h'^''^' dt' 
(4.17) V Jq 

< CKVke^"^ 

and the lemma is proved. □ 
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5 The linear equation 

We continue the analysis of the linear equation (j4.ip . As seen in (j4.5p , when 
expanded in Fourier series, this equation couples the coefficients of indices 
p and —p. The case of indices +1 and —1 is studied in the previous section. 
Using the notations 

(5.1) V = vxe'^ + v-xe~"^ + v" 

we consider the equation ()4.ip for functions with vanishing Fourier coeffi- 
cients of indices ±1 : 

(5.2) L'lU" = F", 

which reduces to the analysis of equations (14. 5p for Fourier p ^ 1. 
The symbol of Lk{dt,p) is 

(-T + mp - k'^p^ -Eq \ 

-T + mp + A;V 'Eq 

k^p^Eo k'^p^Eo k^p^-ry 

which is of course equal to the symbol (j2.3p with with ^ = kp, = —mp, up 
to a change of sign in the first line. 

Assume first that p > 1. In this case, we consider Lk{dt,p) as a pertur- 
bation of 

(idt + mp — k'^p^X 
idt + mp + k'^p'^ 
d! + A:V / 

For the wave operator, we use the classical estimates: 

Lemma 5.1. There is C > 0, such that for all k > 1 andp > \, the solution 
n of 

(5.5) d'^n + k'^p'^n = g, n(0) = dtn{0) = 
satisfies 

(5.6) kp\n{t)\ + \dtn{t)\ < C||5|Ui([o,t]). 

For the Schrodinger equations, we use the following estimates. 
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Lemma 5.2. There are C > and ko > 1, such that for all k > ko and 

p >2, the solutions of 

(5.7) (idt + mp± k'^p^)e = f, e(0) 
satisfy 

(5.8) kV\e{t)\ + \dte{t)\ < C{\\f\\mio,t]) + Wdtflk^im) + l/WI) 
Proof. Standard energy estimates imply that 

\e{t)\<C{\em + \\f\\LHm))- 
Differentiating in time the equation, we obtain 

\dte{t)\ < C{\dtem + \\dtf\y(^0A))- 
The initial condition in ()5.7p implies that 9je(0) = — z/(0). Therefore, 

\ikV ± mp)e{t)\ + \dte{t)\ < C(||/||z.i([o,t]) + II^JIIlHIo,*]) + 1/(0)1 + \f{t)\) 



Recall that m is linked to k through ()3.3p . Thus mp < k'^p + kp and 
f^2p2 _ > ^2 |-p2 _ _ > cfc^p^ for all p > 2 if /c is large enough. □ 

Proposition 5.3. Consider the equation (j4.5p ifii/i initial data 

(5.9) ep(0) = ep(0) = np(0) = 9*^^(0) = 

Then, for p > 2, k > ko, there holds for t £ [0, 1].' 

A;V|ep(i),ep(t)| + \dtep{t) , dtep{t)\ + kp\np{t)\ + \dtnp{t)\ 
(5-10) <C(ll/p, /pllLi([o,t]) + W^tfp, dtfp\\Lmo,t]) 

+ \fp{o)Jp{o)\ + \fp{t)Jp{t)\ + Uphnm))- 

Proof. The lemmas above imply that the left hand side is estimated by the 
right hand side plus 

C(\np{t)\ + ||np, dtfip, k^p^ep{t), A;VepllLi([o,t])) 

The first term is absorbed in the left hand side by kp\np{t)\ for k large 
enough. With Gronwall's lemma, this implies ()5.10p for t G [0,1], with a 
larger constant C. □ 
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When p = 0, there holds: 
Lemma 5.4. When go = 0, the solution of (j4.4p with vanishing initial data 



IS 



(5.11) no = 0, eo(t)= / h{t')dt' . 

Jo 

With the estimates ()5.10p . one deduces the following result 

Corollary 5.5. There are ko and C such that for all k > kg, K , T > 0, 

and all {f",g") with go = 0, satisfying for t G [0, T] 

k"2\\f"{t)\\Hs + k~'2\\dJ"{t)\\Hs < Ke''^', 
\\9"{t)\\H^ < KA;3/4e2'^*. 

the solution of (j5.2p with vanishing initial data satisfies 

k\\e"{t)\\Hs+2 + k--^\\dte"{t)\\Hs < CKe^''\ 
k^\\n"{t)\\Hs + k-^\\dtn"{t)\\Hs < CKe'^^K 

Proof. By Lemma 15.41 there holds 

(5.12) k\eo{t)\ + \dteo{t)\ < Ci^i^e^-* + Vke^^'-^'^e^"'' dt'^ < CKe 
Next, Proposition 15.31 implies that e" satisfies 



2at 



k^dy{t)\\H^ + \\dte"{t)\\H^ 



(5.13) 



< Ck((1 + e^-*) + 1^ A:3/4ga(t-t')e2at'^^/^ < k'/^CKe^'^'. 



Together with (j5.12p this implies the first estimate. 
Moreover, Proposition 15.31 implies that n" satisfies 



k\\n"{t)\\Hs + \\dtn"{t)\\Hs 

Ck((1 + e^^^*) + 1^ fc3/4^.(t-t')g2.t'^^/^ < k'/^CKe 

Since no = 0, this implies the second estimate . . 



(5.14) ^ 



□ 
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6 End of proofs 

First, we note that Proposition 13.31 is an immediate consequence of Corol- 
laries 14.41 and 15. 5i 

It remains to prove Proposition 13. 4i With U = (e,n) and U* = (e*,n*), 
there holds 

(6.1) Afk{U,U*) = if,g)with 



(6.2) f = ne* + n*e, 

(6.3) g = k^dl{Re{ee*)}. 

Proposition 13.41 follows from the next estimates. 

Lemma 6.1. There is a constant C, independent ofk, such that 

(6.4) Vk\\f{t)\\H^ + -^mmH^ < Ce''''\\U\\EHT)\\U*hHT), 

(6.5) ||ff(i)lk^ <Cfc3/V^i[/|| 
Moreover, the mean value go of g vanishes. 

Proof. The first estimate follows directly from the definitions and the in- 
equality 

WabWH" < C\\a\\H''\\b\\H''- 
Next, we note that for e = eie'^ -|- e' and e* = e|e*^ -|- e*' 

dliee*) = dl(e'e*')+Tidlie*'e-'') + et5|(eV^). 

Hence, in norms, there holds 

\\dl{ee*)\\H^ < Wdle'Wns {\\e*'\\ + ||5,e*'f) + ||9|e*'||H» (||e'|| + \\dee'f) 
+|ei|(||5,V'|| + ||e*'||)+|et|(||9|e'|| + He'll) 

and (16. 5p follows. 

In addition, the 0-mean value go vanishes since 5 is a ^-derivative. □ 
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